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1 Introduction 

In this paper we establish the long-time behaviour of the solutions to the discrete three dimen- 
sional Schrodinger and Klein-Gordon equations. We extend a general strategy introduced 
by Vainberg [16], Jensen-Kato [7] and Murata [12] which concerns the wave, Klein-Gordon 
and Schrodinger equations, to the discrete case. Namely, we establish the smoothness in the 
continuous spectrum and the Puiseux expansion for the resolvent of the generator. Then the 
long-time asymptotics can be obtained by means of the (inverse) Fourier-Laplace transform. 

We restrict ourselves to a "nonsingular case", in the sense of [12], where the truncated 
resolvent is bounded at the ends of the continuous spectrum. This holds for a generic 
potential and allows us to get the decay of order ~ which is desirable for applications 
to scattering problems. 

First we consider the 3D discrete version of the Schrodinger equation. 



iip{x, t) = Hip{x, t) := (-A + V{x)) ^'(x, t) 



x G Z^ teR, (1.1) 



^Supported partly by FWF grant P19138-N13, DFG grant 436 RUS 113/929/0-1 and RFBR grant 07- 
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where A stands for the difference Laplacian in Z^, defined by 

\y-x\=l 

for functions ip : —>■ C 

Definition 1.1. Denote by V the set of real valued functions on the lattice 7? with finite 
supports. 

For the potential V, we assume that V The Fourier-Laplace transform 

oo 

i!{x,uj)= I e''^*tlj{x,t) dt, lmcu>0 (1.2) 



reduces f 1 1.11) to corresponding stationary equation 

{H - = -iipo, lmu;>0. (1.3) 

Therefore, 

where R{uj) = {H — uj)~^ is the resolvent of the Schrodinger operator H. 

We are going to use the functional spaces which are the discrete version of the Agmon 
spaces p. These spaces are the weighted Hilbert spaces = with the norm 

= 11(1 + x2)'^/2m||^2 , (reR. 

Let us denote 

B{a, a') = C{ll, ll,), B{a, a') = Oil © ^ © ll) 

the spaces of bounded linear operators from l\ to l\i and from /^©/^ to l\i respectively. 
We write K = Op{K{x,y)) for the operator with kernel K{x,y), i.e., 

{Ku){x) = J2Kix,y)u{y), x G 

Note that the continuous spectrum of the operator H coincides with the interval [0, 12], 
and the resolvent has singularities at points Uk = 4:k, A; = 0, 1, 2, 3 (see [4J). Our main results 
are as follows. For "a generic potential" V &V (see Definition 13.31) . we obtain the Puiseux 
expansion for the resolvent at the singular spectral points Uk- 

R{uJk + uj) = Rko + 0{V^), uj^O, (1.5) 

in the norm B{a; —a) with a > 7/2. Then for initial data ipQ G with a > 7/2 we obtain 
the following long-time asymptotics: 

n 

^/'(■,t) = Y,Cje-''^'^Uj + t/o(t)0± + r±(t), t ±oo. (1.6) 
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Here Uj are the eigenfunctions of operator H, corresponding to the discrete eigenvalues /ij, 
Uo{t) -dynamical group of free Schrodinger equation, 0± G - asymptotic scattering state, 
and ||r±(t)||,. =0(|t|-V2). 

For the proof, we develop the arguments similar to [7], [I2]. Lemma 10.2 of Jensen-Kato 
[7j plays the crucial role in verifying the decay (11. 6p 

We obtain the asymptotics (11. 5p from the invertibility of operator 1 + AkoV, where Ako 
are the coefficients with in the expansion of free resolvent (see. (12.121) . (12.201) ). Note 
that operator 1 + A^qV is invertible for "a generic potential". Asymptotics (11.51) implies the 
boundedness of the truncated resolvent at the singular points of the continuous spectrum. In 
the other words, this points are neither eigenvalue nor resonance for the operator H. Using 
the methods of [12] one can easy proved that the boundedness of the truncated resolvent is 
equivalent to invertibility of the operator 1 + A^qV. 

We also obtain similar results for the discrete Klein-Gordon equation: 



-ipix, t) = (A - - V{x)) ipix, t) 



X G Z^ t eR. (1.7) 



Set *(t) = {^{■,t),ip{-,t)), *o = (^o,7ro). Then ([TTD takes the form 

i^{t) = H*(t), t e M; *(0) = *o, (1-^ 

where 

The resolvent R(co') = (H — cj)^^ of the operator H can be expressed in terms of the resolvent 
R{w), and this expression yields the corresponding properties of R(a;). In particular, we 
obtain the asymptotic expansion of type (11.51) for R(u;), and also the long-time asymptotics 
of type (II. 6p for the solution of (11.80 . 

Let us comment on previous results in this direction. Eskina [3J, and Shaban and Vain- 
berg [H] considered the difference Schrodinger equation in dimension n > 1. They proved 
the limiting absorption principle for matrix elements of the resolvent and applied it to the 
Sommerfeld radiation condition. However, [3l [H] do not concern the asymptotic expansion 
of R{uj) and the long-time asymptotics of type (II. 6p in the operator norms. 

The asymptotic expansion of the matrix element of the resolvent R{u) at the singular 
points Uk was obtained by Islami and Vainberg pj. They used this expansion to prove 
the long time asymptotics for the solutions of the Cauchy problem for the difference wave 
equation. The main feature which differs the present paper from [4j is that here all asymptotic 
expansions hold in the weighted functional spaces, not on compacts as in [Ij. In fact, the 
asymptotic expansion of the resolvent (11.50 in the norm of B{a, —a) was the main technical 
challenge in this paper. Additional difference is that the long time asymptotics here is 
obtained for the Schrodinger and Klein-Gordon equations, contrary to the wave equation in 

The asymptotic expansion of the resolvent and the long time asymptotics (11.60 for hy- 
perbohc PDEs in W (continuous case) were obtained earlier in pjj, ^15j, [16j, [17j, and 
for the Schrodinger equation in [7], [H], [S], [12]; also see [T3] for an up-to-date review and 
many references concerning dispersive properties of solutions to the continuous Schrodinger 
equation in various norms. 
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The results of present paper extend the results of [9] and [TOj from difference ID and 2D 
equations to difference 3D equations. 

The paper is organised as follows. In Section [2] we derive the asymptotic expansion of the 
free resolvent. The limiting absorption principle and the Puiseux expansion of the perturbed 
resolvent is proved in Section [HI In Section H] we prove the long-time asymptotics fll.6p . In 
Section [6] we extend the results to the discrete Klein-Gordon equation. 

2 The free resolvent 

We start with an investigation of the unperturbed problem for the equation f 1 1.11) with 
V{x) = 0. The discrete Fourier transform of u : — > C is defined by the formula 



u 



After taking the Fourier transform, the operator Hq = —A becomes the operator of multi- 

3 3 ^ 

plication by 0(6*) := 6 — 2 ^ cos 9j = 4 ^ sin^ y: 

i=i i=i 

- ^{6) = <j){e)u{0), e e (2.1) 

Thus, the operator Hq is selfadjoint and its spectrum coincides with the range of the function 
0, that is Speci^o = S := [0, 12]. Denote by Ro{uj) = {Hq—lu)^^ the resolvent of the difference 
Laplacian. Then the kernel of the resolvent Ro{uj) reads 

JioK = 5^ /S^'*«.-eC\E. (2.2) 

r3 

Lemma 2.1. The free resolvent Rq{u!) is an analytic function o/cj G C \ S with the values 
in B{a, a') for any a, o' E M. 

Proof. For a fixed G C \ S, we have (j){6) — a; 7^ for G T^. Therefore, (f){6 + iC,) — 00 
for 6 G T^, C, G M^, if ^ 7^ is sufficiently small. Hence, the function l/{(j){6) — uj) admits 
analytic continuation into a complex neighbourhood of the torus of type {6 + iC, : 6 E T^, ^ E 
: 1^1 < 6(00)} with an 6{uj) > 0. Therefore the Paley- Wiener arguments imply that 

Ro{u;,x-y) < C((5)e-^l^-^l 

for any S < S{uj). Hence, Ro{uj) is the Hilbert- Schmidt operator in the space B{a, a'). □ 

2.1 Limiting absorption principle 

Now we are interested in the traces of the analytic function Ro{uj) at the cut S. We can 
write 

Roiuj ± te, z) = F-^^-^^ = — / -— de, z G Z3 (2.3) 







with e > 0. Note that the hmitinff distribution — — is well defined if u is not a 

critical value of the function 6(6), i.e. the level line 6(6) = uo does not contain the critical 



points with V(j){6) = 0. The critical points 6 = {6i, 62, 63) can be easily calculated: 6i = /cjvr, 
ki G Z. Therefore, the critical values are u = 0,4,8, 12. For the free resolvent the following 
limiting absorption principle holds 

Proposition 2.2. For a > 3/2 the following limits exist as e 0+.' 

B(o-,-<t) 

RQ{uj±ie) ^ i?o(cu±iO), G S\ {0,4,8,12}. (2.4) 

Proof. First we prove the convergence (12.41) for any fixed z, follow [3]. Let Xj{^)y 3 = 1, ^ 
are the sufficient small partition of unity on the torus T^. Then 

± ^) = E 5^ / = E ± (2.5) 

where Dj = supp Xj- If {0(^) = u;} fl = 0, then the function Rqj{uj ± ie, z) is continuous 
for £ > and 

|i?0j(^±«e,2)| < Q < 00, zGZ^ £ > 0. (2.6) 
Now let Sj = {(j){6) =uj}n Dj. Then G Dj 

6 = s + tn(^s), 

where s G S'j, and n{s) is the external normal vector to Sj at the point s of unit length. Let 
us introduce the new variables {s,t). Then 

/ , • N 1 /" -is... /" Xj{s + tn{s))e-''''^'^'J{s,t) 

Roi(uj±te,z) = — - e ''""ds / — — ^ " , dt. (2.7) 

'■'^ ' ' J J t^{s + tn{s)) T ie 

Sj -a{s) 

where J{s,t) is the Jacobian, ip is the smooth function, and a(s), b{s) > 0. Note that 
J(s,t)|(=o = 1; and ■iIj{s + tn(^s))\t=Q = |V0(s)|. We will prove the following lemma: 

Lemma 2.3. Let (p{t, z) is the smooth function satisfies 

\^{t,z)\<C, \dt^{t,z)\<C\zl tG [-5,(5], ;2GZ3, (2.8) 

and let ip^t) is the smooth function such that iplt) 7^ z/t G [—5,5]. Consider F{±e,z) = 
5 

Vit, z) 

— — dt. Then 

tip{t) =F ie 

F{±e,z) ^ F{±0,z), ,£^0+, ^GZ^ (2.9) 



and 

ee(0,l] 



sup |F(±£,2)| < C(ln(l + 1^1) + 1), zeZ^. (2.10) 
e(o,il ^ ^ 
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Proof. Let us rewrite F{±e, z) as 

5 5 5 

F(±e, z) = (^ 0, z) / -— (^ 0, z) / , ^ ,;\ 7 + / \i . dt. 

~S -5 -5 

Then, 

Fi±e,z)^Fi±0,z)^±^n^-^ [ t^^l^clt+ [ ^^^^l^^dt (2.11) 



^(0) m J tm J tm 

-5 -5 

as £ — >■ 0+. By ( 12.81) the first and the second summand in RHS of (12. lip can be estimated 
by the constant which does not depend on \z\ G Z^. Let us estimate the third summand 
in RHS of (12. lip . For \z\ < 1/6 this summand also can be estimated by the constant. For 
\z\ > 1/(5 we obtain by (12.81) that 

s 

1 



/l^^^l- / - / 



< --C\z\+C\n\z\ < C\n\z\ 

\z\ 



-5 \t\<l/\z\ l/k|<|t|<5 

Lemma is proved. □ 
Lemma 12.31 imphes that 

loj{uJ ±xb,z) — ^ nojl 

and 



Roiicu ± ie, z) Roj(uj ± iO, z), £ — > 0+, zEZi^, 



sup \Roj{uj ±ie,z)\ < Cj( \n{l + \z\) + 1 

6(0,11 ^ 



£6(0,1] 

Evidently, that the all resolvent Rq satisfies the similar properties. Hence by the Lebesgue 
dominated convergence theorem 

^ (1 + \x\'^)-''\Ro{uj±ie,x -y) - Rf){uj ±iO,x - y)\'^{l + \y\'^)~'' ^ 0, e 0+ 

with a > 3/2. Then the Hilbert-Schmidt norm of the difference Ro{u ± ie) — Rq{uj ± iO) 
converges to zero. Proposition 12.21 is proved. □ 

Remark 2.4. Differentiating (\2.3\) with respect to u, we obtain similarly that for G S \ 
{0, 4, 8, 12} the derivatives dj^Ro^u ± iO) belongs to B{a, —a) with a > 3/2 + k. 

2.2 Asymptotics near singular points 

We obtain the asymptotics Ro{^^) near the singular points Uk. We consider "eUiptic" points 
uJi = 0, UJ4 = 12 and "hyperbolic" points = 4, c^a = 8 separately. 



2.2.1 Asymptotics near elliptic points 

Here we construct the Puiseux expansion of the free resolvent Rq{uj) near the point cui = 
(the expansion near the point = 12 can be construct similarly). 

Proposition 2.5. Let N = 0, 1,2... and a > N + 3/2. Then the following expansion holds 
in B{a, —cr): 

N 

R^(uj) = ^AifcCU^/2^0(cj(^+i)/2^, \uj\ ^0, argcuG (0,27r). (2.12) 

k=0 

Here Aik E B{a, -cr) with a > k + 1/2. 

Proof. The resolvent -Ro('^ =t iO) is represented by the integral (12.21) . Fix < 5 < 1 and 
consider < |cj| < 5^/2. We identify with the cube [— 7r,7r]'^ and represent Rq{uj,z), 
z = X — y, a.s the sum 



57r3 / 6(9) — uo Stt^ / <Z)iy)— u; 



2^3 -pa 
i?0l(c^,2;) + -Ro2(^,^), 



where x(^) is the smooth function and 

X{r) = 

Since (f){9) = |6'p + (9(|6'|'^), then Ro2{^, z) is analytic function of ui in |a;| < 5^/2, and 



0, r > 2(j 

1, < r < a 



|5ii?02(^,^)|<^^^^, kl<5V2, ^GZl (2.14) 

Hence it suffices to prove the asymptotics of type fl2.12p for _Roi- Let us choose the system 
of coordinate in which ^3 coincides with the axe z and rewrite Rqi{uj, z) as 

25 

n..(-.^) = ^ [ U^^^^f^W (2.15) 

8vr^ y (j)(rn) - uo ' 

\n\=\ 

Here r = |^'|, Q' = rn, 4>{6') = (t>{d). Note that is even function. Hence, 

2(5 2(5 

/2oi(-,.) = ^ / f 7 ^(,)e-Nn3,2 ^ ? x(yN-r^cir x^ ^ ^^^^^^ 

25 

x(r)e~*'''^'"3r^(ir^ 



-25 



{rn) — uj 



ds. 
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where 5"+ = {|w| = 1, n^, > 0. (We extent x(^) on M as even function). Let us apply the 
Cauchy residue theorem to the inner integral in fl2.16p : 



Roi{uj,z) 



p{uj,n) 



ds + 



St 



S+ [-2S,S]U[5,25] 



x{r)e ^''l^l^sr^ dr 
f>{rn) — uo 



ds (2.17) 



+ 




S+ Tg 



^~ir\z\n3^2 

(j){rn) — u) 



ds = -Roi(w, z) + Rliiuj, z) + RliiuJ-, z) 



Here q = q{uj,n) is the solution to (f){qn) = in the low half-plane Img < 0, p{uj,n) = 
dr4>{rn) , = {|r| = 5, Imr < 0}. Since (f){rn) = + ^{rn) with $(m) = O(r^), then 



r=q 



q{uj, n) = -^/iu{l + Q{uj, n)), p{q, n) = ^/uJ{2 + P{uj, n)), 



(2.18) 



where Q{u!,n), \l/(co',n) are the analytic function of u, and Q{uj,n), "${00,11) = 0{uj). Let us 
consider the third summand in the RHS of (ETTj): 



2n vr/2 



Rl^{uj,z) 



in-" 



dp J da 

Ts 



(2.19) 



2tt 



tt/2 



in^lzl 



dp J dr 





irdr 
(r,a,/3) - iu 



d{e 



-ir\z\ cos 



Since |0i(r, a,/?) — > C{5)5, then the operator value function Rq^{ijj) is analytic and all 

its derivatives in respect to u are bounded in i3(cr, —a) with a > 1/2 for \u!\ < 6^/2. Hence 

i?Q^ admits an expansion of type 02.121) . 

The expansion of -Rqi can been obtained similarly. 

Finally, let us consider the first summand in RHS of (ETTj): 



27T tt/2 



■^TT^lzl 



dp 





2n 



q{uJ,7r/2,P) 

p{uj, n){q sin a + daq cos a) Stt^I^I L7 p{uj,7i/2,P) 



dp 



27r it/2 

- ! dp I e-^^l"l™^" d- 



2-K 



q sm a 



l + Q(cu,7r/2,/3) 

p(a;,n)(gsina + c^agcoso;)] Stt^I^;! L7 2 + P(a;, 7r/2, /3) 
00 

2-n 7r/2 

+ dp g-iv^(l+Q('»^>a,/3)k|cosa ^ 



dp 



{l + Q{u;,a,P)f 







(2 +P{uj, a, p)){l + Q(w, a, p) +daQiuj, a, /?)/ tana) ' 



Here q{uj,n) = q{uj,a,P). Expand the integrand, we obtain the asymptotics of type (12.121) 
for -Ro2(i^, z), since Q, P ^ uj. □ 

Remark 2.6. The expansion (12.121) can been differentiated N + 1 times in B{a, —a) with 
a > N + 3/2: 

N 

dlRoioj) = ( ^ ^fccu'^/^j ^ C'(^(A^+i)/2-r-)^ 1 < r < AT + 1. 



A;=0 



2.2.2 Asymptotics near hyperbolic point 

Here we construct the Puiseux expansion of the free resolvent Ro{uj) near the point to'2 = 4 
(the expansion near the point = 8 can be construct similarly). The main contribution into 
(12.21) is given by the corresponding critical points (0, 0, vr), (0, tt, 0) and (vr, 0, 0) of hyperbolic 
type. 

Proposition 2.7. In B{cr, —a) with a > 2N + 7/2 the following expansion holds: 

N N 

i?o(4 + cj) =^A2feCj'= + V^^52fctu'^ + C(cj^+^), |cu| ^ 0, Imcj > 0, (2.20) 

fc=0 k=0 

where A2k,B2k G i3(cr, — a), wits a > 2k + 3/2 are the operators with kernels A2k{x — y), 
B2k{x - y)- 

Proof. For c<j = 0^2 = 4 the denominator of the integral 02.21) vanishes along the curve 
0(6*) = 4. We will study main contribution of points (0,0, vr), (0,7r,0) and (tt, 0,0) of the 
curve which are critical points of (j){0). The contribution of other points of the curve can 
be proved by methods of Section 12.11 For example, let us consider the integral over a 
neighbourhood of the point (7r,0,0). Let (,{9) be a smooth cutoff function, equal 1 in a 
neighbourhood of the point (7r,0,0) (the other properties of C,{9) we specified below). For 
Imcu > denote 

^, , If e-^("i''i+^2e2+^3e3) 
Q{uj,z)- 



8773 J 4){e) ~A-uj 

<^iei+^2e2+z,e,)^I^Q-^ d^i (102 dOs e-'^^^ r e-i{zie{+z2e2+z,,es)^^i^0,-^ ^q,^ 



Stt^J 4sin2| + 4sin2|-4cos2|-cu Sn^ J 4 sin2 | + 4 sin^ | - 4 sin^ f - a; ' 

where 9[ = Oi-n, 9' = (6*;, 6*2, 6*3), and Ci{9') = ({9). We suppose that Ci{9') is symmetric 
in 9[, 6*2 and 6*3. Then the exponent e~*'^^i^i+^2^'2+^3^'3) ^g^j^ be substituted by its even part, so 
we have 



00 00 00 



, _ e— f f fcosiz,9[) cos(^2^2) cos(z3^3) Cii0')d0'id02d9s _ , , 

J J J 4sin2% +4sin2% -4sin2^ tt^ 

000 ^ ^ ^ 

(2.21) 

Let us obtain the expansion of type (12.201) for Qi. We change the variables: si = 2sin y, 

S2 = 2sin^, and S3 = 2sin ^ and choose the cutoff function ( such that Ci(^') = C2(|'5p), 
with smooth function (2- Then 



00 00 00 




000 



F{ z,sj,sl,sl)C2{\s\'^) dsi ds2 ds3 

-|- ^2 — "S^ — U! 



Qi{uj, z) 
where 

F{z, s\, si, S3) = cos(zi5((si)) cos{z2g{s2)) cos{z3g{s3))J{sl)J{sl)J{sl). 
We change the variables again: 



si = r cos ip,S2 = T sin (p, S3 = S3. 



Then 



— si — UJ 








where 

7r/2 



j F{z, cos^ ip, sin^ yj, Sg) di{) 





We change the variables once more: 

p^ = T^ -s\ = cos 2ip, p2 = 2rs3 = sin 2ip, (2.23) 

where R, ip are the polar coordinates on the plane (r, S3). Then |pp = Pi + P2 = -^^5 hence, 
IpI = R'^, and 

r= = ^. .1 = ^. (2.24) 

Therefore, dpidp2 = 4|p|c?rc?S3, and (12.221) implies 

00 

Q,{u;,z) = j [j (^|^^|^'j|^| v^^^^^Pi)^P2^ (2-25) 

R 

where /i(|p|,pi, 2;) = ^1(2;, ^^^^^^ , ^^^^^^ )^2(|p|)/4-\/2. Now we can specify all needed properties 
of cutoff function: 

suppC2(|p|) n {p G : P2 > 0} C n = {(pi,p2) : -5 < Pi < 5, < p2 < 5} 

with some < 5 < 1. We consider < |ct;| < 5/2, Imcj > 0. Denote r = r(p) := |p|. The 
function /i(r, pi, 2;) can be expanded into the following finite Taylor series with respect to pi: 

h{r,pi,z) = ho{r,z) + hi{r,z)pi + ... + hN{r,z)p'^ + H^ir, pi, z)p^ , (2.26) 

where hk{r, z) are polynomial in z of order 2k, and 

\HM{r,pi,z)\<C\z\^'', \d,,HM{r,pi,z)\<C\z\^^^\ (pi,P2) G [-5,5] x [0, 5], r = |p|. 

(2.27) 

^iep Let us substitute (I2.26P into (12.251) and consider the summands with hk{r,z), k = 
0,1, ...,iV: 

hkj r, z)p\^r^ pi fhk{r,z)^r + pi f fc-i , ^^ , 

7 dpidp2 = / Pi H-cupi +...+UJ H ]dpidp2 

[pi-uj)r J r \ pi-ujJ 

n n 

fc-i 



X: a.,(.)^^- Mllilv^r^rf^^rf^^ (2.28) 

f-^ J [pi- uj)r 



yakj{z)uj^ + UJ / ^dr / ^ — , 

^—'^ n n 



ii 



where akj{z) are polynomial in z of order 2k. Let us calculate the integral 
Vl + cost/' dtp } 2V2c/(sin|) f y^dt 



cosV'-cu/r J 1 - 2 sin^ I - cu/r 7 - ^ 





log ^^^ + ^^=, lmu;>0. (2.29) 



Here y/r > 0. Function z = i/r — is analytic in Imcu > with the values in Imz < 
0, Re2; > 0, and function ( = logw is analytic in |w| < 1, Imw > 0, where log(— 1) = vri. 
Substitute (12:291) into (ESHD, we get 



/■ hk{r,z)p\ ^4 t ^ j , k f I ■ . ^ V'^^r'' \ hkir,z)dr 

J {p -iu)r ^^ Pidpidp2 = ak,j{z)uj-^ +uj J \m-\og ' 



n J=o \ ^ + V~/ 

(2.30) 

Let us expand /iA;(r, 2;) into the following finite Taylor series with respect to r: 



hkir, z) = hk,oiz) + h,iiz)r + ... + h^N-kizy'-^ + Hk,N-k{r, zy'-\ (2.31) 

where hkj{z) are polynomial in z of order 2{k + j), and \Hk,N~kir, z)\ < C\z\'^'^ , <r <6. 
The following lemma is true 

Lemma 2.8. Let < \uj\ < 5/2, Imcu > 0. Then 



m - log ^ = siicu) + Gcu'v^, (2.32) 

where si are analytic in < \uj\ < S/2, Imu > 0, Ci E M. 

We shall prove this lemma in Appendix B. Now fl2.30p . (12.311) and (I2.32p imply that for 
< \uj\ < 5/2, Imcu > 



/h ir z") ^ iv N-k 

/ ^ ^/r + pidpidp2 = ^dkj{z)uj^ + lo^^ ^ bk,j{z)uj^ + d^+i^^ 

n ^^'^ i=o i=o 



j=0 j=0 

^ ' 1 / r—uj \ , ^ nT 1 



h II . V 2r \ Hh Ar_u(r, z)r dr 
+ 00^ \m-\og ^-j= ''^ , (2.33) 



\ ' V 2r 



where \dk,,{z)\ < C{k,j)\z\^^, \bk,j{z)\ < C(A;, j)kl'('+^'^ and {d^+iA^, z)\ < C(iV + 
l,k)\z 



2N 



Let us consider the integral in RHS of f l2.33p for < leu | < 5/2: 

2\uj\ 5 



_ logilii (2.34) 



In the first summand we change the variable: r = Ic^It. Then 

1 - 



















'\uj\ 


/( 










^ 



iri — log ■ 



1 + 



t—uj/\lu\ 
27 



T — UO \Ul!\ 



(2.35) 



2Af|, .\N-k 



Let us consider the second summand in RHS of (12 .34^ : 



N-k 



,N-k^ 



Hk,N-kir,z)r^ ^ ^/2(^do + c?i^ + -- + c^7V-fc^^ + (^"jv-fc(c^/r)^^) 



dr 



2|a;| 
<5 



Hi 



'k,N-k{r, z) + d.ujr^'-^-^l'' + ... + rf^v-fc^^-'r-^/^) dr + Un-kiyJ. z) 



2|^| 



MAr_fc(CJ,Zj 





N-k 

Ujiz)uJ^ + MAr-fc(a;, z), 

j=0 



(2.36) 



2Af|, ,|Af-fc 



where |uj(2;)| < C|z| |; |uAr_fc(u;, z)!, \uN-ki^, z)\ < C\z\ \uj 
Step ii) It remains to consider the contribution into (12.251) from the remainder H]^[r, pi, z)pi : 

CHN{r,pi, z)p^dpidp2 



(pi - uj)r 



+ pi 



HNir,pi,z)^r + pi / 



Pi 



7V-1 



N 



n 

TV-l 

5^^,(z)c.^- + a;^ 

j=0 



pl-UJ 



dpidp2 



HNjr, pi, z)y/r + pi dpidp2 
(pi - uj)r 



(2.37) 



where |wj(z)| < C|2;p^. For the integral in RHS of (12.371) the following estimate is true 
Lemma 2.9. Let < Ic^l < 5/2, Iracu > 0. Then 

HNir, pi, z)y/r + pi dpidp2 



(pi - uj)r 



We shall prove this lemma in Appendix C. 
Step Hi) Finally, we obtain 



< C\z\''' \n' \z\, \z\ > 1. 



(2.38) 



N 



N 



Qi{uj,z) = ^qk{z)u'' + y/uj^pk{z)u'' + qN{uj,z), \uj\ 0, 



k=0 

\2N i„2 I ^1 1 IN 



k=0 



where |g}v(cu,z)| < C\z\^^ \n' \z\\u\^ . Further, pk{z) = 0{\zr), and qk{z) = 0{\z\^^) for 
< k < N. Therefore, qk{z) = (9(|-2p^), since qk{z) do not depend on A^. □ 
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Remark 2.10. The expansions (12.201) can be differentiated 2N + 2 times in B{a, —a) with 
a > 2N + 7/2. More precisely, 

N N 

d:Ro{4 + uj)=d:(^Yl + E + (^(^""^''l , 1 < r < 2iV + 2 

k=0 k=0 

3 Perturbed resolvent 

3.1 The limiting absorption principle 

In the next proposition we develop the results of [3], [E] for the 3D case and prove the 
limiting absorption principle in the sense of the operator convergence. It will be needed for 
the proof of the long-time asymptotics fll.6l) . 

Proposition 3.1. Let V E V , a > 3/2. Then the following limits exist as e — >• 0+ 

R{uj±ie) ^^^^ R{uj±iO), G S \ {0, 4, 8, 12}. (3.1) 

Proof. The existence of the limit (13. ip follow from the representation R = Ro{I + V^-Ro)"^, 
Proposition 12.21 and from the invertibility of the operators I + VRq{uj ± iO. The invertibility 
of the operators I + VRq{uj ± iO follow from [HI Theorem 9]. □ 

Remark 3.2. For u; G S \ {0,4,8,12} the derivatives d^R{ijj ± zO) belong B{cr,—a) with 
a > 3/2 + A;. 

3.2 Asymptotics near singular points 

In this sections we are going to obtain an asymptotic expansion for the perturbed resolvent 
R{uj±iO) near u = 0,4,8,12. 

Definition 3.3. i) A set W gV is called generic, if for each V &V we have aV G W, with 
the possible exception of a discrete set o/a G M. 

a) We say that a property holds for a "generic" V , if it holds for all V from a generic subset 
ofV. 

Theorem 3.4. Let a > 3/2. Then for "generic" V & V the following expansion holds: 

R{uj) = Dl + 0{^/^), ^ 0, argcj G (0,27r) (3.2) 

in the norm of B{a, — o"), where Di is the operator with the kernel Di{x,y). 
Proof We use the relation 

R{iu) = T-\uj)Roiuj), where T{uj) := I + Ro{uj)V. (3.3) 
According to (l2rT2il . 

T{lj) = I + AioV + 0(y/II^), \uj\^0, argcj G (0,27r). (3.4) 



J.4 



Let us prove that for "generic" V G V the operator T{uj) is invertible in for sufficient 
small \uj\ > 0. It is suffices to prove that for "generic" V eV the operator T(0) = Op 6{x — 

y) + AioViy) is invertible in Z^^, or the operator 



Op (1 + x'r'^/\6{x -y)+ A,,V{y)){l + y') 



is invertible in Let us consider the operator 



A{a) = Op[(l + a;2)--/2 [S{x ~ y) + aAioV{y)) (1 + y^y/^] = 1 + a/C, a G C. 



2W/2] 



For a > 3/2 



y) = (1 + x^)-'^/Miol^(2/)(l + y'r^' G /^(Z2 X 



Hence, K{x,y) is a Hilbert-Schmidt kernel, and accordingly the operator }C = Op{K{x,y)): 
P is compact. Further, A{a) is analytic in a G C, and ^(0) is invertible. It follows 

that A{a) is invertible for all a G M outside a discrete set; see |2j. Thus we could replace the 
original potential V by aV with a arbitrarily close to 1, if necessary, to have Tr{0) invertible. 
Now (13. 3p and fl3.4l) imply that for sufficiently small > 



R{u;) = {I + T(0) + 0{V^)y\A,o + 0{V^)) = T{0)-'A^o + O(v^). 



(3.5) 



□ 



Remark 3.5. i) The expansion of resolvent near the second elliptic point uj 4 = 12 is similar 
to the expansion ( (3.21) . 

a) The expansion of type ( (3.2j) near the hyperbolic points uj2 = 4: and = 8 require larger 
value of a. Namely, for "generic" V eV 



R{uJk+uj) =Dk + 0{V^), 



\UJ\ 



0, ImLc;>0 



in the norm of B{(T, —a) with a > 7/2. 

Hi) These expansion can be differentiated two times in B{a, —a) with cr > 5/2 for elliptic 
points and with a > 7/2 for hyperbolic points. In these cases d'^R{tOk + to) = O^cu^^^"^), 
k = 1,2,3,4. 



4 Long-time asymptotics 

Theorem 4.1. Let a > 7/2. Then for "generic" V eV the following asymptotics hold 



0{t 



-3/2n 



, t — > 00. 



(4.1) 



B{u-a) 



Here Pj denote the projections on the eigenspaces corresponding to the eigenvalues fij G 
M\[0,12], j = l,...,n. 
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Proof. The estimate (14. ip is based on the formula 

^-UH ^ — L I e-'*'^R(uj)duj, C > m8ix{12]\fiA, j = l,...,n}. (4.2) 
2m J 

Therefore 



n 

-itH 



■'^^ [0,12] [0,12] 

The main contribution into the long-time asympotics gives the integrals over the neighbour- 
hoods of the singular points = 0, 4, 8, 12. For example, let us consider the integral over 
the neighbourhood of the point uoi = 0. The expansion fl3.2p and Remark 13.51 imply 

pW(cj) = 0(0'' V^), ^ +0, G M, = 0, 1, 2 (4.3) 

in the norm of B{a, —a) with a > 7/2. Then by Lemma 10.2 from ^ we get 



in the norm of B{a, —a) with cr > 7/2. Where C{uj) is the smooth function, supp ( G (—1, 1), 
and C(^) = 1 for ^ ^ [0, 1/2]. The integrals over the neighbourhoods of the points a; = 4, 8, 12 
can be estimated similarly. □ 



5 Application to the asymptotic completeness 

We apply the obtained results to prove the asymptotic completeness which follows by stan- 
dard Cook's argument. 

Theorem 5.1. i) Let tpo G Then for "generic" V E V for solution to (11. 11) the following 
long time asymptotics hold 

n 

^(■'^) = EE^'^.^"'*''''"^. + U,m±+r±{t), t ^ ±oo, (5.4) 
i=i kj 

where u^^ are the corresponding eig en functions, (p± G are the scattering states, and 

\\r±{t)\\i2 ^ 0, t^ioo. 

a) Furthermore, 

||r±(t)||,. = 0(|trV2) (5.5) 

ifipo G ll with a > 7/2. 



Proof. We will consider the asymptotics for t ^ oo. It suffices to prove the asymptotics 
( 15.41) for ipQ G with a > 7/2 since is dense in and the group Uq is unitary in P. Let 
us apply the projector P onto the continuous spectrum to both sides of fll.ll) : 

iP^j = HoPiJ + VP%Ij. (5.6) 

Hence, the Duhamel representation gives, 

t 

Pip{t) = t/o(t)P^(0) + j Uo{t - T)VPip{T)dT, t G M. (5.7) 



We can rewrite (15.70 as 

oo oo 

P^(t) =f/o(t)(p^(0) + J Uo{-r)VP^{r)dT'^ - J Uo{t - r)VP^{r)dr = + r+{t). 

t 

(5.8) 

oo 

Now we can prove that the ffist integral in (15.81) converges, and the function J Uq{—t)V Pip {T)dT 



belongs to P. Indeed, the unitarity of Uo(t) in implies by ( 14. ip that 



f/o(-r)yP^(r)||prfr = J ||FP^(r) ||,2rfr < C J ||P^(r) ||,2_^dr 



oo 

< c f r^/^umwpdT <c 







Here we used that P is bounded operator in since P = I — ^ Pj- The estimate ( 15.51) for 
r+{t) follows similarly. □ 



6 The Klein-Gordon equation 



Now we extend the results of Sections [3M] to the case of the Klein-Gordon equation (11.81) . 
The resolvent R(ti;) can be expressed in term of the resolvent R{uj) as 



R(cu) 



-i(l + cu2p(cj2-m2)) cjP(cu2-m2) 



The representation ( 16. 9p implies the following long time asymptotics: 
Let a > 7/2. Then for "generic" 1/ G V 



± i=i 



—itv- 



0{t 



-3/2^ 



OO. 



B((T,-cr) 



(6.9) 



Here are the projections onto the eigenspaces corresponding to the eigenvalues v- 
±^m2 + cjj, j = 1, . . . ,n. 



± 



7 Appendix A 

Let us prove Lemma 12. 81 For / = we get 



.5 



^ " ' ~ \ dr ^ , / . . ^~ \ ~ 



h = I I vrz - log ^--= = 2Vr - uj m - lo. 



2V2u: t , = 2VS^ (m - log ' \K ] - .2V^log f " 



\ V 2<5 



^ — l-iJf , — ^ 
= so{u) — i2\'2uj log — ^ — TTv^u = so{uj) + CoytJ, 

where Sq, Sq are the analytic functions of to, Cq = —n\/2. Further, for / > 1 we get 



^ f 1 — ^ / ^-T^ \ rr-lrlrr. / 1 



. V \ rdr I r- / . ^ Y 25 

TTi — log —j= —j=^ = /O V — TTZ — log , 



1^ / -| /t W \ 1^ 

- 211 ; ^ «-log dr + 2v^w / 



\ ' Y 2r / 

.5 



-2//i + 2/cj/j^i + 2y2cj [ ^(r^-^ -ujr^-^ + ... + (-ujy-^ +^^-^) 

J y/r\ r + UJ/ 



si{uj) - 2lli + 2luli-i + Si(cu) - i2V2^{-ujf log ^ - 2'kV2^{-u)\ 

1 ~|~ 2 ^ 



5 

where 5;, 5; are the analytic functions of u. Hence, Ii = si{uj) + CiuJ'-Jijj, where si is the 
analytic functions of lo and Ci G M. 

8 Appendix B 

Here we prove Lemma I2.91 We estimate only integral over H+ = {0 < pi,p2 ^ i^}- The 
integral over H \ H_|_ can been estimated similarly. Let us split the integral over H_|_ into two 
integrals: 



{pi-iu)r J {pi-uj)r 

n+ 

f HN{r, \uj\,z)y/r + pidpidp2 _ j , j 

+ / 7 ^ — -'1 + ^2- 

(Pi - uj)r 



n 

Similar fl2.29p we obtain 

6 



J I I ■ 1 V 2'' \ HN{r, \uj\,z)dr 

^2 = / K« - log 



UJ 



J.5 



Note that 



log 



2r — \Jr — uj 



log 



r + 



2r + -Jr — UJ 

< I log |r + ct;| I + 2| log \ + \/r — uj\\ < 2| log \r — \uj\\\. 
Then ([2IH]) imphes 

I J,| < / l±}}^it^dr < C\zr. 



UJ\ 



Further, for \z\ > 1 let us split Ji as 



Ji — Jii + Ji2 + J: 
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where Jn is integral over Hi = {{pi, P2) G n+ : |r| < J12 is integral over 112 = 

{(pi,P2) e n+ \ni : |pi - < and J13 is integral over lis = n+ \ (Hi un2) (see 

picture 1). 



v = -p>|co|/2 




|co|/2 l®l |co|+v 
Figure 1: \lu\ - l/\z\^/^ < 0. 
By (12.271) and the inequality \pi — u!\ > \pi — \u!\\ we get 

< C\z 



Pi 



I Jill < 

Hi 

7t/2 1/|2|*/^ 



Ipi - (^\r 



\pi - \uJ\\^yr + pidpidp2 ^ ^|^|27v+2 f + Pidpidp2 



Hi 



7r/2 1/|2|*/^ 7r/2 l/l^l*/^ 







For the second integral we obtain similarly 



|27v+2 y y/r + pidpidp2 ^ (^|^|2JV+2 y 



|2Ar 



n2 n2 
since l/^/r < \z\^/^ for (pi,p2) e Hs, and iHa] < 25/\z\^/^. Finally, (E^ZD implies 



|Ji3|<C^kP^ / '^Pl^^^<C\zr\n'\z\ 



n 



since for any vertical interval / G Ha we get 

/ -7==^ = HP2 + J Pi + Pi) <Cln \z\. 



References 



[1] Agmon S.: Spectral properties of Schrodinger operator and scattering theory, 
Ann. Scuola Norm. Sup. Pisa, Ser.IV 2, 151-218 (1975) 

[2] Bleher P.m., On operators depending meromorphically on a parameter, Moscow 
Univ. Math. Bull. 24, 21-26 (1972) 

[3] EsKlNA M.S., The scattering problem for partial-difference equations, in Mathematical 
Physics Naukova Dumka, Kiev, 1967 (in Russian), 248-273 

[4] ISLAMI H., Vainberg B., Large time behavior of solutions to difference wave opera- 
tors, J. Commun. Partial Differ. Equations 31, no. 1-3, 397-416 (2006) 

[5] Jensen A., Spectral properties of Schrodinger operators and time-decay of the wave 
function. Results in L'^{W^), m > 5, Duke Math. J. 47, 57-80 (1980) 

[6] Jensen A., Spectral properties of Schrodinger operators and time-decay of the wave 
function. Resuhs in ^^(R^), J. Math. Anal. Appl 101, 491-513 (1984) 

[7] Jensen A., Kato T., Spectral properties of Schrodinger operators and time-decay of 
the wave functions, Duke Math. J. 46, 583-611 (1979) 

[8] Jensen A., Nenciu G., A unified approach to resolvent expansions at thresh- 
o\ds,Reviews in Math. Physics 13, no. 6, 717-754 (2001) 

[9] KOMECH A., KOPYLOVA E., Kunze M., Dispersive estimates for ID discrete 
Schrodinger and Klein-Gordon equations, J. Appl. Anal. 85, no. 12, 1487-1508 (2006) 

[10] KoMECH A., KoPYLOVA E., Vainberg B., Dispersive estimates for 2D discrete 
Schrodinger and Klein-Gordon equations, J. Funct. Anal. 254, no. 8, 2227-2254 (2008) 

[11] Lax p., Phillips R., Scattering Theory. Academic Press, New York (1989) 

[12] MuRATA M., Asymptotic expansions in time for solutions of Schrodinger- type equa- 
tions, J. Funct. Anal. 49, 10-56 (1982) 



[13] ScHLAG W., Dispersive estimates for Schrodinger operators: A survey, preprint 
,math.AP/0501037, 

[14] Shaban W., Vainberg B., Radiation conditions for the difference Schrodinger oper- 
ators, J. Appl. Anal. 80, no.3-4, 525-556 (2001) 

[15] Vainberg B.: Behaviour for large time of solutions of the Klein-Gordon equation. 
Trans. Moscow Math. Soc. 30, 139-158 (1974) 

[16] Vainberg B.: On the short wave asymptotic behaviour of solutions of stationary prob- 
lems and the asymptotic behaviour as t oo of solutions of non-stationary problems, 
Russ. Math. Surveys 30, 1-58 (1975) 

[17] Vainberg B.: Asymptotic Methods in Equations of Mathematical Physics. Gordon 
and Breach, New York (1989) 



